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Abstract 

The standard derivation of Schrodinger's equation from a Lorentz- 
invariant Feynman path integral consists in taking first the limit of 
infinite speed of light and then the limit of short time slice. In this or- 
der of limits the light cone of the path integral disappears, giving rise 
to an instantaneous spread of the wave function to the entire space. 
We ascribe the failure of Schrodinger's equation to retain the light 
cone of the path integral to the very nature of the limiting process: it 
is a regular expansion of a singular approximation problem, because 
the boundary conditions of the path integral on the light cone are lost 
in this limit. We propose a distinguished limit, which produces an 
intermediate model between non-relativistic and relativistic quantum 
mechanics: it produces Schrodinger's equation and preserves the zero 
boundary conditions on and outside the original light cone of the path 
integral. These boundary conditions relieve the Schrodinger equation 
of several annoying, seemingly unrelated unphysical artifacts, includ- 
ing non-analytic wave functions, spontaneous appearance of disconti- 
nuities, non-existence of moments when the initial wave function has 
a jump discontinuity (e.g., a collapsed wave function after a measure- 
ment), the EPR paradox, and so on. The practical implications of the 
present formulation are yet to be seen. 

Keywords: Lorentz-invariant path integral, subluminous propagation, bound- 
ary conditions 
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1 Introduction 

The original derivation of Schrodinger's equation from Feynman's path in- 
tegral postulates a classical non-Lorentz-invariant action pQ, [2], [3], [I]. 
The Lorentz-invariant (discrete) Feynman path integral [5], which propa- 
gates within the light cone emanating from the support of the initial wave 
function, is usually approximated by expanding the Lorentz-invariant action 
in powers of x/c first (c is the speed of light) and then expanding the path 
integral in powers of the time slice At. In this approximation the light-cone 
disappears and Schrodinger's equation in the entire space is recovered, as in 
the non-relativistic case (see details below). If the order of limits is reversed, 



the discrete path integral converges to the initial wave function and does not 
propagate. This anomaly in the approximation of the path integral can be 
ascribed to the observation that the aforementioned power series is a regular 
expansion of a two-parameter singular perturbation problem. Therefore, the 
construction of a uniform asymptotic approximation to the Lorentz-invariant 
Feynman path integral calls for a singular perturbation approach. 

In this paper we construct a uniform asymptotic approximation to the 
Lorentz-invariant Feynman path integral by identifying a distinguished limit 
of large c and small time slice At, such that cAt — > 0, but cvAt — > oo. More 
specifically, as the speed of light c is constant, and explicitly present in the 
problem, the time slice At cannot be assumed arbitrarily small, as explained 
above. In fact, the distinguished limit requires that At be bounded below 

by At 3> — - ph 6.2 x 10~ 22 sec (for an electron). In this distinguished 

mc 2 
limit the leading order approximation to the path integral is a solution of the 

initial value problem for Schrodinger's equation with zero boundary condi- 
tions on and outside the light-cone emanating from the support of the initial 
wave function (see Figure [1]). This subluminous Schrodinger equation can be 
considered an intermediate model between non-relativistic and relativistic 
quantum mechanics and it is valid for the above mentioned time resolution. 
We stress, however, that it is not our aim to construct a comprehensive 
relativistic quantum mechanics theory (in, say, the sense of [6]). 

Notably, the emergence of the light cone relieves Schrodinger's equation 
of several annoying, seemingly unrelated mathematical and unphysical arti- 
facts. These include non-analytic wave functions, spontaneous appearance of 
discontinuities [7], non-existence of moments when the initial wave function 
has a jump discontinuity (e.g., a collapsed wave function after a measure- 
ment) [8], the EPR paradox [9] (relativistic quantum mechanics theory not 
withstanding), and so on. On the other hand, conventional quantum mechan- 
ics is recovered inside the light cone after an appropriate relativistic delay. 
The practical implications of all this are yet to be seen. 



2 Superluminous propagation 

The Feynman path integral formulation of quantum mechanics in one dimen- 
sion postulates that the propagation of a free particle is defined by the action 



functional 

t 



S(x(t))= I ^^-dt (2.1) 



o 
as a convolution with the propagator 

exp|^(x(t))}. (2.2) 

This means that the propagation of the time-sliced wave function at time 
t + At is given by [I] 

oo 

. s / tti f f im (x — y) 2 1 . , . 

^ t + At) = \/^iKAtJ ^{it-im}^^ (2 - 3) 

— oo 

The limit At — > is found PQ, [1] by changing the variable of integration to 
y = x — z^/fiAt/m, which converts ( 12. 3ft to 



V>(M + At)=M/i /exp{^}^(x-W^,t) <fe (2.4) 



-oo 
oo 



2«/ exp T^ x 



iz 2 



dz 



,hAt , . , 1 2 frA£ , . , 

ip(x, t) - z\ ip x (x, t) + -z ipxx{x, t) + 

m 2 m 

= ^{x, t) - - ^ xx (x, t) + 0[[ 

2 m \\mj 

hence Schrodinger's equation follows. 

If instead of the non-relativistic Lagrangian C(z) = z 2 /2 in ( 12. ip . we 
use the unique relativistic Lagrangian of a free particle in one dimension 
C(z) = — \J\ — z 2 [5], it leads to the Lorentz invariant action 

t 

S(x(t)) = I ' mc 2 C f^lj dt'. (2.5) 

o 



Now, the Lorentz invariant propagation is given by 

ip{x,t + At) =M~ l x (2.6) 

mc 2 At 



exp < i 



V c 2 At 2 



n 

where the normalization constant M is given by 



H[c 2 At 2 -(x-y) 2 ]iP(y 7 t)dy, 



oo 

M = cAt f exp {^^C (^f\ } H [c 2 At 2 - (x - yf] dy. (2.7) 



Note that modifying C(z) to C(z) = 1 — yi — z 2 does not change the integral 
equation (j2.6p . Obviously, the wave function ip(x, t) defined by (12. 6p vanishes 
on and outside the light cone. 

Taking the first the limit c — > oo converts ( 12. 6 ft to ( 12.31) . and then the limit 
At — > 0, as in ( 12 ,4p . recovers the Schrodinger equation in the entire space 
and the zero boundary conditions on light cone disappear and so does the 
light cone. Evidently, the wave function (12. 3p spreads to the entire line in- 
stantaneously, thus rendering Schrodinger's equation superluminous. Taking 
the limits in the reverse order does not lead to a wave equation. 

The decay of the integral in ( 12.31) for large \x\, similarly to the Fourier 
transform, is determined by the regularity of ip(y, t). Thus jump discontinu- 
ities in ip(y,t) lead to a decay of ip(x,t + At) as (x)" 1 , which prevents the 
existence of moments, energy, and other artifacts [10] . 

3 Subluminous propagation 

From the mathematical point of view, the approximation of the integral (1 2 . 6 [) 
by taking limits in the order described in Section [2] is a regular expansion 
of a singular approximation problem. The singularity of this expansion is 
manifested in the loss of the boundary conditions on the light cone. We 
propose here a distinguished limit that removes the singularity and the the 
ensuing artifacts. Specifically, we consider a general Lorentz-invariant action 
functional (1 2 . 5 j) . where the Lagrangian C(z) is an analytic function near the 
origin with 

£(0) = 0, £'(0) = 0, £"(0) = 1. (3.1) 



The assumption that the action (12.51) is Lorentz-invariant implies that the 
time-sliced wave function it defines cannot propagate outside the light cone 
emanating from the support of the initial wave function. We postulate, as 
above, the integral equation 

ip{x,t + At) =M~ l x (3.2) 

oo 

/ exp { j ^ir lc (jm) } H ^ At2 - {x ~ y) ^ ^ t] dy > 

— oo 

with the normalization factor M given by 

oo 

A/" = cAtje W !^^c(^^H[c 2 At 2 -(x-yr}dy 

— oo 

1 
= cAt/exp^M^, (3.3) 

-1 

where 

« = =^. (M 

To investigate the behavior of ip(x,t) on the boundary of the light cone 
emanating from the initial support, we rewrite (13. 2p in the explicit form 

x+cAt 

^x,t + At)=Af- 1 J exp|^^/:r^)U(y,*)dy- (3-5) 

x—cAt 

If Xq is a point on the boundary of the initial support, we assume the light 
cone emanating from it at time At is x = Xq — cAt. At this point (13.51) is 

if>(x - cAt, At) = A/" -1 x 

XO 

J ex P {^£( I °- c 'f;-^ }^.0)^ = 0, (3.6) 



x -2cAt 
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Figure 1: The thick line is the initial wave function, normalized in L 2 (M). The 
light cone for c = 1 and t < 2 is the domain above the x-axis and above the black 
region in the (x, t) plane. 



because the entire (open) interval of integration, (x — 2cAt,a;o), is outside 
the support of ip(y, 0). Note that ip(xo — cAt, At) = 0, even if ip(xo, 0) ^ 0. 
This boundary condition persists as long as the line x = Xq — ct is on the 
light cone (see Figure [1]). Thus, if (a, b) is a finite interval on the x-axis 
outside the support of ip(x, 0), then the wave function ip(x, t) vanishes in the 
triangle in the (x,t) plane, whose vertices are the points [a, 0], [6,0], and 
[(a + 6)/2, (6 — a) /2c]. This implies, for example, that if in the two slits 
experiment (a, 6) is the interval separating the two slits, there will be no 
interference pattern between the slits before time (6 — a) /2c. Furthermore, 
the influence on the wave function of every finite interval in the support of 
the initial wave function is confined to the light cone emanating from the 
interval (see Figure [Q . 

It is easily seen that the solution tp(x,t) of fl3.5[) is Lorentz-invariant on 
the discrete time grid and it is supported on the light cone emanating from 
the initial support. 



3.1 The distinguished limit 

Equations (I3.3P and ( 13. 4 p suggest the distinguished limit of large c and small 
At such that cAt — > 0, but cyAt —¥ oo (in appropriate dimensionless units). 
In more precise mathematical terminology the meaning of this distinguished 
limit is the assumption that £ — > oo, or more specifically, 

h . L 2 m , _. 

< At < — — , (3.7) 



mc 2 Ti 

where L is a characteristic length, e.g., the width of the initial support of the 
wave function (the lower bound is due to the presence of c, and the corre- 
sponding characteristic time and the upper is needed for the approximation 
of a finite difference quotient by a derivative, see below). 

For a point (x, t) inside the light cone emanating from the initial support, 
we change the variable of integration to y = x — zcAt and write (13. 2p as 

?/>(£, t + At) = A/^cAt x (3.8) 



/ exp < iC(z) — - — > ip(x — cAt(l + z),t) dz. 
-l 

If the light cone containing (x, t) is the interval Xq — ct < x < Xq + L + ct, we 
assume that x — ct + e = x < x + L + ct for some positive e independent of 
At. Then integration in (13.81) extends over the interval — 1 < z < —1+e/cAt, 
which contains the entire interval — 1 < z < 1 if cAt is sufficiently small. Now 
we expand in Taylor's series to convert (I3.8P to 

i 

/( mc 2 At 1 

expUC(z) — - — ^x (3.9) 



-i 

1 



ip(x,t) — zcAtip x (x,t) H — z 2 c 2 At 2 ip xx (x 1 1) + 



dz. 



The left hand side of the inequality ( 13. 7ft ensures that £ ^> 1 (see ( I3.4p ). so 
the integrals in (13. 9p can be evaluated by the stationary phase method [TT] . 
The right hand side of ( 13. 7p serves to ensure that ip(x,t + At) — ip^x^t) = 
Atipt{x,t) + o(At). Thus the expansion (13. 9 j) with the assumptions ( 13. ip . 
( 13. 7p gives 

# t (x, t) = -^Pxx + O (r 1/2 ) , (3.10) 



which in the (formal) limit £ — > oo is Schrodinger's equation in every interior 
point of the light cone. 

Because the support of ip(x, t + At) cannot exceed the support of ip(x, t) 
by more than cAt on either side, Schrodinger's equation A3. 10|) with a finitely 
supported initial wave function has to be solved inside the light cone ema- 
nating from the initial support (see Figure [TJ with vanishing boundary con- 
ditions. In particular, zero boundary conditions have to be imposed on the 
triangles mentioned above. Thus, for example, if the initial wave function 
vanishes on a finite number of finite intervals and outside the minimal finite 
interval containing its support, the zero boundary conditions on the triangles 
are given only for a finite time, but the zero boundary condition on the light 
cone emanating from the minimal interval containing the support are given 
for all times. The Schrodinger equation has to be solved piecemeal outside 
the triangles for finite time intervals and then in the entire light cone em- 
anating from the above mentioned minimal interval. Obviously, the initial 
conditions at times corresponding to every apex of a triangle, as mentioned 
above, is the wave function obtained from solving the Schrodinger equation 
in the light cones emanating from the support up to that time. Thus, after 
the time corresponding to the highest triangle the Schrodinger equation has 
to be solved in the light cone emanating from the minimal interval, with the 
initial value of the wave function constructed up to that time. An important 
result of this structure is that at all times the wave function is analytic inside 
its (finite) support and has finite moments at all times. 

We apply the above procedure to the case discussed in Section [2] C(z) = 
1 — y/1 — z 2 . In this case all integrals can be expressed explicitly in terms of 
Bessel and Struve functions [12] without invoking the saddle point expansion. 
The Lorentz invariant propagation is now given by (12. 6p . so the change the 
variable of integration to y = x — zcAt and expansion in Taylor's series 
converts (12. 6 j) to 

i 
ip(x,t + At) = M^cAt I exp | -i Vl - z 2 - 1 m ° - \ x (3.11) 



-i 

1 



ip(x,t) — zcAttp x (x,t) + -z 2 c 2 At 2 ip xx (x,t) + ■ 



dz. 



The coefficient of — ip xx (x,t) in (13. lip is shown in Figures |2] and |3j In 
m 




Figure 2: The lower curve is the real part and the upper is the imaginary part of 
the coefficient of £rip xx (x,t) in (15TT1) for < £ < 20 
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o- VfcA^AAAAAA/VVVUVV\/\AA.VVAAAWWIAAA 



Figure 3: The lower curve is the real part, oscillating around its limit and upper 
curve is the imaginary part, oscillating arou: 
of %ip xx (x,t) in QHZETD for 1000 < f < 1200. 



curve is the imaginary part, oscillating around the limit value of the coefficient 5 
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the limit At — > 0, that is, £ — >■ 0, all coefficients in the expansion ( 13. lip 
(after dividing by At) converge to 0, giving ip t (x,t) = so the path integral 
converges to the initial wave function. However, the limit cAt — > and 
£ — > oo in the expansion (13.111) yields (13.101) . It is evident that the path 
integral converges to the solution of the subluminous Schrodinger equation 
as At —?• 0, except for a narrow range of At, where the restriction (13. 7p is 
violated. 

4 Quantization 

In the presence of a potential $(x) the action (12.51) is replaced with 



t 



2rl X ( t ' 



S(x(t)) 



and the exponent in (13. 2p becomes 



1 I _ 2 r ( X ~ V 



mc 2 C -^^- - $(x{t')) 



dt', (4.1) 



s r £ U^~* (!/) J Ai ' (42) 

The Schrodinger approximation is now 

h 2 

m t {x, t) = -—ijj xx (x, t) + q>(x)if>{x, t)+0 (r 1/2 ) (4.3) 

within the appropriate light cones confines. Quantized energy levels appear 
when the potential has a proper local minimum. To see how discrete energies 
emerge in the subluminous Schrodinger equation, we consider the elementary 
case of the harmonic oscillator with potential $(x) = ^mu 2 x 2 and assume 
that the support of ip(x, 0) = <f(x) is the interval (— x ,x ) and that the 
boundary conditions are ^>(±(xo + ct),t) = 0. First, we nondimensionalize 
equation f!4.3[) by setting 



X = V\I > i = -> i>(x,t) = U(y,r), e = -J—, (4.4) 

moo oo c V m 
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which converts the initial and boundary value problem for Schrodinger's 
equation (14. 3p into 

2iU T (y, r) = -U w (y, r) + y 2 U(y, r) for \y\ <Vo + ^ (4.5) 



U(y,0) = v L^lLy 1/(±( WJ + I),r)=0. (4.6) 

Introducing the fast and slow time scales a = r/e and r, respectively, and 
setting U(y, r) = V(y, a, r), we rewrite (14. 5 j) as 

2i(V T + e" 1 ^) = V ro + 3/V (4.7) 

and assume the regular outer expansion [TT] 

y = V^ + eV 1 + --- . (4.8) 

The resulting hierarchy of equations 

V a ° = 0, 2iV2 = 2iV T -V y ° y + y 2 V ,... (4.9) 

implies that l/ is independent of e and a and thus cannot satisfy the 
boundary condition on the light cone y = ±(yo+cr). The solvability condition 
for (14.91) . obtained from averaging with respect to the fast variable a, is the 
Schrodinger's equation 

2iV T -V y ° y +y 2 V° = 0, (4.10) 

which has to be solved on the entire line with the initial condition (j4.6p . The 
first approximation V° is thus merely the scaled classical quantum harmonic 
oscillator 

oo 

V°(y,r) = V°0M) = ^a„e-*'/ s i(x), (4.11) 

n=0 

where E n are the energies, if) n (x) are the (real valued) eigenfunctions of the 
harmonic oscillator on the entire line and a n are the coefficients in the ex- 
pansion of the initial wave function V°(y,0) = <p(x) in the eigenfunctions 
ip n (x). 

12 



Higher order terms in the outer expansion (14. 8 p cannot satisfy the bound- 
ary conditions (14. 6p . so a boundary layer correction [IT] is needed at y = 
±(yo + t/e). We thus introduce the boundary layer variable and function 
near the right boundary 



T 

yo + --y 

e 



U(y,r) = W( V ,r) 



(4.12) 



2iW° 



-W° nv + t 2 W° 



Schrodinger's equation (14.31) gives 

2ie 2 W T + 2iW v = -W m + (ey + r - e 2 7]) 2 W. (4.13) 

Expanding W -= W° + eW 1 + ■ • • , we obtain the boundary layer equation 

(4.14) 
and the boundary and matching conditions 

W°(0,t) = -V°(y + ^,r^, W°(oo,r) = 0. (4.15) 

The solution is given by 

W°(t], t) = -V° (y + -, t) exp |-r/ i + Vr 2 - 1 j. (4.16) 

With a similar construction at the left end of the interval, we get the uniform 
expansion for e< 1 



U(y,r) 



H 



T 



(y -yo) - -n 



V°(y,r) 



(4.17) 



-V°(yo + -,r)exp 



+ V~i 



-V° i^-y ,rjexp 



+ Vi 



t + e(y - y) 



r + e(y + y) 



Equation (I4.17P indicates that, for large r/e, the expansion of the wave 
function in the eigenmodes of Schrodinger's equation on the line is recovered 
inside the light cone, away from the boundary. To regain a single mode, the 
initial wave function has to be of the boundary layer form ( 14.171) at time 
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r = and V°(y, 0) has to be the an eigenfunction of Schrodinger's equation 
on the entire line. 

Obviously, energy is preserved, because boundary conditions and integra- 
tion by parts of the energy integral give 

XQ + Ct 

~ •■ -icipipt{xo + ct) + icipipt{-xo - ct) - i / [^) t ^t + i>i>u] dx = 0. 



dt 

— XQ—Ct 

The energy at time t is up to errors of order of magnitude of the contribution 
of the boundary layer 



xn+ct 

oo 






, m,n=0 

Xq—CZ 



-a^th^e-^h^^MMx) 



dx 



XQ+Ct 



J2 a m a n E m e i{En - Em)t/n / i> m {x)Mx) dx 

^,n=0 •* , 

' —XQ—Ct 

OO 

-* 2_^i \ a ™\ 2 Em as ct —)• oo. (4-18) 



5 Discussion 

The standard non-relativistic Schrodinger model of quantum mechanics con- 
sists of the initial value problem for Schrodinger's equation in the entire 
space. In the conventional interpretation of this model any boundary condi- 
tions imposed on the wave function can be represented in terms of appropriate 
potentials in Schrodinger's equation. In this paper we adopt an alternative 
model, which consists of a Lorentz-invariant Feynman path integral that per- 
force vanishes on and outside the light cone emanating from the support of 
the initial wave function. This boundary condition does not require any po- 
tentials, rather, it is an inseparable part of this model of quantum mechanics. 
As explained in this paper, the path integral converges to a limit as the time 
slice At — y 0, except in a narrow range where the restriction (13.71) is violated. 
The distinguished limit cAt — y and £ — y oo of the path integral is the 
solution of the subluminous Schrodinger equation (the solution that satisfies 
zero boundary conditions on and outside the original light cone of the path 
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integral). We can view, therefore, the subluminous Schrodinger equation 
as an intermediate model between non-relativistic and relativistic quantum 
mechanics. The results of the standard non-relativistic quantum theory are 
recovered in the intermediate model, albeit after a relativistic delay. 

To clear the issue of how the results of this paper relate to existing theo- 
ries, we point that 

1. This paper mainly addresses issues related to the original (non-relativistic) 
Schrodinger equation. 

2. Ours is not a relativistic quantum theory, but rather a conventional 
one that upholds the basic tenets of physics: all propagation must be 
confined within a relevant light-cone. 

3. Why should one care about such matters as light-cones in the non- 
relativistic regime? The quick answer is that neglecting the speed of 
light and the resulting boundary conditions on light-cones is, mathe- 
matically speaking, a singular perturbation, which results in unphys- 
ical artifacts such as the disappearance of moments and energy when 
the initial wave function is discontinuous (e.g., after a collapse due to 
a measurement) or the instantaneous spread of information into the 
entire space. These artifacts merely reflect the mathematical short- 
comings ingrained in the original Schrodinger formulation. 

4. Recall that similar difficulties occur in diverse physical theories. For 
instance, air viscosity may be practically negligible, but if it is neglected 
altogether, gas dynamics predicts no drag force on airplanes wings. For 
no matter how irrelevant it is elsewhere, close to a moving body there 
is a boundary layer, where viscous effects are dominant. Thus, while 
viscosity manifests itself locally, its effect is global! 

5. A closer kin to our problem is the heat equation and its unphysical 
artifact of instantaneous spreading of initially localized perturbation to 
the entire line. This artifact can be attributed to the disregard of the 
acoustic speed limitation in the derivation of the diffusion equation. 
If the acoustic bound is enforced on the process, finite propagation 
fronts emerge and the unphysical artifact of instantaneous propagation 
disappears [13J. 
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6. Perhaps the most succinct summary of our discussion would be to say 
that it is not enough that we know the limitation of our model; the 
model itself should contain its own limitations, as is the case at hand: 
the Schrodinger equation with boundary conditions contains the limi- 
tation on propagation. 

7. The non-relativistic propagator of Schrodinger's equation has some of 
the properties of the Fourier transform: the decay at infinity of the wave 
function reflects the smoothness of its initial value. Thus a jump discon- 
tinuity is converted by the propagator into a decay of the wave function 
at infinity as O da;!" 1 ), which although square integrable, has no mo- 
ments [10]. In contrast, none of this anomaly exists for Schrodinger's 
equation with zero boundary conditions on and outside the light cone. 

To conclude, while confining the Schrodinger equation within a relevant 
light cone may only mildly extend the scope of its applicability, which is yet 
to be seen, it definitely relieves it of annoying paradoxes and artifacts. Given 
the central position of Schrodinger's equation in modern physics, our exten- 
sion should be expected to have implications beyond mere scientific esthetics. 
The basic predicament of relativity manifests itself in our theory only in cru- 
cial 'junctions' of space-time: the theory limits the domain of propagation, 
endows the wave function with a definite front, thus eliminating the unlim- 
ited precursors and causing a time delay for the conventional quantum effects 
to re-emerge. 
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